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We  study  the  limiting  distribution  of  the  amoimt  of  charge  left  in  some  set  by  an 
infinite  ^stem  of  charged  Markovian  particles,  when  the  charge  distribution  belongs 
to  the  domain  of  attraction  of  a  symmetric  a-stable  law.  The  limits  are  symmetric 
a-stable  generalized  random  fields.  Thdr  multiple  int^rals  are  built  in  a  similar 
manner.  We  also  study  the  renormalizability  of  these  families  of  random  fields  and 
use  the  construction  to  simulate  stable  fields  on  and 


1.  INTRODUCTION 


Stable  processes  have  been  studied  intensively  in  recent  years.  Samorodnitsky 
and  Taqqu  (1990),  Weron  (1984),  Rosinski  and  Woyczynski  (198j),  Rosinski  (1986), 
Kallenberg  and  Szulga  (1989),  ozulga  (1992),  Maejima  (1990),  Janicki  and  Weron 
(1991),  are  some  of  the  raerences  related  to  this  work.  Because  all  these  authors  work 
on  an  abstract  measiurable  $pace,  or  with  specific  modeb  on  R},  we  were  motivated  to 
define  families  of  stable  random  fields  on  R^,  d  >  1,  and  of  generalized  stable  fields  on 
the  Sobolev  space  We  bmld  the  fields  as  limits  of  sums  of  functionals  on  paths  of 
Markov  processes  and  show  that  using  this  construction  one  can  easily  analyze  some 
of  the  properties  of  the  fields  in  the  l^t. 

Motivated  by  work  of  Maejima  (1990)  on  self-similar  stable  processes,  we  choose 
to  study  the  renormalizability  of  the  families  of  stable  fields  and  their  functionals. 
The  property  of  renormalizability  of  families  of  random  fidds  was  defined  in  Adler 
and  Epstein  (=Feldman)  (1987)  and  Epstein  (1989),  as  an  extension  of  the  property 
of  self-similarity  to  families  of  fields.  These  authors  discuss  the  renormalizability  of 
fEunihes  of  Gaussian  fields.  Here  we  show  that  stable  fields  built  &om  self-similar 
Markov  processes  are  renormalizable. 

Our  method  of  construction  of  the  stable  fields  also  allows  us  to  simulate  these 
fields  on  a  computer  and  display  the  surfaces  to  see  what  they  look  like.  Our  moti- 
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vation  in  this  direction  comes  from  work  of  Janicki  and  Weron  (1991)  on  simulation 
of  stable  processes  (all  on 

We  now  look  at  our  results  from  a  somewhat  different  point  of  view.  In  recent 
years,  much  attention  has  been  given  to  the  description  of  infinite  systems  of  particles 
moving  according  to  some  law  (usually  Markovian).  Some  of  these  papers  e.g.,  Adler 
and  Epstein  (1987)  and  Adler  (1989),  deal  with  particle  systems  whicli  behave  as 
follows:  Initially  (at  time  zero)  a  number  of  independent  particles  pop  into  existence 

at  locations  within  the  space  according  to  a  Poisson  point  proce^  with  intensity 
A.  (Actually  the  initial  distribution  of  the  particles  was  created  differently  in  the 
abow  papers,  but  one  could  use  Pmsson  point  processes  instead.)  The  partides  then 
move  acoarding  to  some  Markov  law.  A  podtive  or'negative  cnarge  is  assigned  to 
each  particle  initially,  according  to  a  Rademacher  random  variable.  The  chtfge  of 
each  decays  exponentially  with  time.  Hie  number  of  particles  in  the  system  is  then 
set  to  infinity  and  the  limiting  distribution  the  charge  left  by  the  system  in  a 
set  in  A^,  after  all  the  partides  have  lost  thdr  charge,  is  studied.  The  limiting  fidd, 
which  is  indexed  by  sets,  or  more  generally  by  functions,  was  shown  to  hsLve  a  Gaussian 
distribution  in  Adler  and  Epstein  (1987).  In  Fddman  and  Eachev  (1993),  the  authors 
obtain  limiting  fields  that  have  sub-Gaussian,  Laplace  and  other  distributions,  by 
changing  the  initial  distribution  of  the  partides  appropriately. 

The  second  aim  of  this  work  is  to  answer  the  question,  “What  happens  to  the 
limiting  charge  distribution  if  the  initial  charge  of  each  partide  follows  a  symmetric 
a-stable  law?” 

As  in  the  above-dted  work,  our  main  tools  are  limit  theorems  for  U-statistics. 
Here  we  use  results  found  in  Szulga  (1992)  on  the  convergence  of  resampled  U-statistics 
to  multiple  stable  integrals. 

This  paper  is  organized  as  follows;  In  Section  2,  we  present  our  construction  of 
generalized  stable  random  fields  and  their  multiple  integr^.  The  proofs  are  given  in 
Section  3.  In  Section  4,  we  show  how  to  construct  stable  fidds  on  A^.  In  Section  5  we 
study  the  renormalizability  of  the  families  of  stable  fidds  and  their  multiple  integrals. 
Section  6  is  devoted  to  simulation  results  for  some  stable  processes  and  fidds  on  A^ 
and  A2. 


2.  STABLE  RANDOM  FIELDS  AND  THEIR  FUNCTIONALS.  g 

< 

Let  us  now  define  the  particle  system  of  the  Introduction  precisely.  On  an  & 

arbitrary  probability  space  (O,  P),  take  n  independent  symmetric  Markov  pro-  g 

cesses  Vi,  V5j,...,Vn  with  values  in  A*^,  each  process  starting  according  to  some  S 

finite  initial  measure  m.  Let  pt(z,y)  be  thdr  common  transition  density  function; 

Pt(®>y)  =  Pt(y»a;),  X,  y  €  A^,  <  >  0  and  /^p«(x,y)dy  =  1,  for  each  i  €  A^.  The 
corresponding  Green’s  function  is  pven  by 

g{x,y)  =  g\x,y),  g^(x,y)=  f  e“^pt(x,y)dt,  (1) 

^0  ns  cRA&i 

riC  TAB 

If  K  is  one  of  the  processes  in  the  sequence,  then  it  can  be  thought  of  as  describing  jannou.  ced 
the  motion  of  a  partide  in  the  system.  We  now  discuss  the  charge  assigned  to  each  istificatior 

of  the  n  particles.  Let  the  probability  space  be  rich  enough  to  support  the  sequences  — - 

of  i.i.d.  random  variables,  y„^i,  Ynfi,  •  •  •  which  are  independent  of  the  V's  and , 
for  which  ’  ’ 
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where  Co  =  cosx)dx/x^'^“.  Let  oi,  02, . . .  ...  be  a  sequence  of  i.i.d. 

Rademacher  variables  (on  the  same  probability  space),  i.e.,  P(<7i  =  1)  =  P{ai  —  ^ 
—1)  =  1/2,  which  are  independent  of  the  V’s  as  well  as  the  V’s.  Assign  to  each 
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particle  V^-  a  charge  We  will  call  o^iVn^i  the  charge  associated  with  the  Markov 

particle  VJ.  The  particle  charges  clearly  belong  to  the  domain  of  attraction  of  a  stable 
law. 

We  now  describe  the  evolution  of  the  system  in  time.  When  particle  i  with  charge 
time  t  passes  throiigh  a  point  i  in  the  space  i2^,  it  leaves  there  a  charge 

Let  A  €  B{R^)  be  a  Borel  set  in  the  space  R^.  We  would  like  to  find  the 
amount  of  charge  left  in  the  set  A  after  all  particles  have  lost  their  charge  and  iii  the 
limit  of  increasing  initial  particle  density,  i.e.,  we  are  interested  in  finding  the  limiting 
distribution,  as  n  — »  oo,  ^  the  normali^  sum 

More  generally,  define  the  bilinear  form 

(/i  h)  =  (/,  f)m  :=  /  m{da)dx  dy  f{x)g\a,  x)g{x,  y)f(y)  (3) 

JR^ 

where  g  and  g^  were  defined  in  (1).  Let  =  Sj{g)  be  the  Sobolev  space  of  C°° 
functions  of  finite  norm,  where  the  inner  pr^uct  in  5^  is  ^ven  by  (3).  For  /  in  Sj 
we  define  the  functional  of  the  space  of  the  paths  of  the  Markov  process  V  by 

Ff(v)  =  r  F-‘f(v(t))dt  (4) 

Then,  we  study  the  weak  convergence,  as  n  -♦  oo,  of  the  finite-dimensional  distribu¬ 
tions  of  the  stun 

Mf)  =  j2<’iyn.iFf{Vi)  (5) 

»=1  t 

Clearly,  if  (|/|,  |/|)  is  finite  for  all  functions  /  from  the  Schwartz  space  of  C°°  functions 
that  decrease  at  infinity  faster  than  any  polynomial,  we  can  take  the  Schwartz  space 
as  the  index  set  of  $n-  This  is  possible,  for  example,  when  the  Markov  process  V  is 
Brownian  motion. 

We  now  define  the  limiting  field.  A  random  variable  X  is  said  to  have  a  symmetric 
a-stable  (SaS)  distribution  if  its  characteristic  function  is  given  by 

Eexp{ttX}  =  exp  (6) 

The  SaS  distribution  is  thtis  characterized  by  index  a,  0  <  a  <  2,  and  the  scaling 
parameter  a.  A  stochastic  process  is  SaS  if  all  its  finite-dimensional  distributions  are 
SaS,  i.e.,  iff  any  linear  combination  of  its  components  is  an  SaS  random  variable.  On 
a  separable  finite  meastue  space  (T,T,ff)  let  X{')  be  a  SaS  Levy  random  measure 
with  control  measure  ft  and  Levy  measure  u(dx)  =  c^^dx/x~^~^y  where  c“  is  as  in 
f2).  This  means  that 

(i)  for  every  A  €  T,X(A)  is  a  SaS  random  variable  with  characteristic  fimction 

Eexp(t<X(A)}  =  exp  {-p(A)|<|®}  , 

(ii)  X(‘)  is  independently  scattered  and  a-additive,  i.e.,  for  disjoint  sets  Ai, . . . ,  A^, . . . 
in  T,  the  random  variables  A’(Ai),...,Jt(Ajt)  are  independent  and  X(U^jAj)  = 

I  (/!;). 


We  will  denote  hy  Xf  the  integral  of  /  with  respect  to  measure  X  and  by  X^  f 
the  multiple  integral  of  order  k  with  respect  to  measure  X.  The  construction  of  these 
is  given  in  Kallenberg  and  Szulga  (1989). 

Now  we  choose  a  specific  separable  finite  measure  space  (T,T,p).  Let  T  = 
D([0,oo))  denote  the  path  space  of  the  ^-valued  Markov  process  V,  let  B{T)  =  T 
dmote  the  Borel  <r-al^bra  of  the  cylinder  sets  in  T,  and  let  p  denote  the  measure 
induced  by  the  process  V.  Thus,  V‘,  Vi,  are  i.i.d.  random  variables  uniformly 

distributed  on  (T,T,p).  Tbke  Jf(*)  to  be  a  SnS  Levy  random  measure  on  (T,T,/j,) 
with  control  measure  ft. 

A  fiunily  d  random  variables  is  called  a  generalized  random  field 

if  the  following  two  conditions  hold: 

(a)  4  is  a  linear  random. fiinctional,  i.e.,  9(af  +  bh)  =  a$(/)  +  6$(h)  a.s.  for  all  /, 
A  €  Sfi,  and  all  a,  6  €  A; 

(b)  $  has  a  version  with  values  in  the  dual  space  Sj. 

Recall  from  (4)  the  definition  of  the  functional  Ff  on  space  T.  We  now  consider  a 
fronily  of  random  variables 


{«(/),  /eSj}S  F,(u)X{du),  f  e  Sj}  (8) 

Proposition  1.  Let  0  <  a  <  2.  Then  (8)  deBnes  the  generalized  stable  random  field 
on  such  that  for  eveiy  /  €  5^  and  t  €  R, 


Eexp{it«(/)}  =exp|-|tr  . 

s«p{-|f|«ElF/(V)|“}  (9) 

Theorem  2.  Let  0<o<2.  Asn-»oo,  the  finite-dimensional  distributions  of 
converge  weakly  to  those  of  the  stable  random  field  {^(/),  /  6 

In  order  to  construct  a  fcth  order  integral  with  respect  to  Levy  random  measure  ^ 
X ,  we  have  to  conrider  systems  of  k  Markov  processes.  For  a  function  /t  of  k  variables, 
define  the  functional 

r.-.r (10) 

JQ  Jo 

Since 

E \Ph  (n.,  ■ .  ■ . Vi.)f  =  (ft,  ft),  i  >  1,  (11) 

where 

{fk^hk)--=  /^^»«(<ia)/t(x)^2(a,x)^(x,y)hfc(y)dxdy  (12) 

P(x,y)  =  ^(a:i,yi)...y(xjt,yi),  y^(a,x)  =y2(ai,zi)...yVt)®ib)  (13) 

(see  Feldman  and  Rachev,  1993),  we  restrict  our  parameter  set  to  fimctions  from  the 
space 

•=  {fk  ■  fk  =  /fc(®l,  •  ••,**)  symmetric  in  C°°(R^*)  with  (|/fc|,  |/jt|)  <  oo} 

(14) 

We  are  interested  in  studying  the  limiting  distribution,  as  n  -»  oo,  of  the  normalized 
sum 

® n(/lr)  •“  •  •  •  ^ii,  •  •  •  Ynit  ^ft  (l^'i  >  •  •  •  i  ) 

i<ii<...<iic<n 


(16) 


Proposition  3.  Let  0  <  a  <  2.  The  Held 

e  sj)  &  h  €  sf} 

. ^ 

is  weU  deBned  and  is  a  generalized  random  £eld  an  Sj. 

Theorem  4.  Let  0  <  a  <  2.  As  n  -*  oo,  the  Bnite-dimensional  distributions  of 
{♦n(/4),  fk  €  5j|  convezige  weakly  to  those  of  ft  €  Sj}. 


3.  PROOFS 

We  will  now  present  the  proofs  of  the  results  stated  in  Section  2.  We  start  with 
the  following  lemma. 

Lemma  5.  For  0  <  a  <  2  and  every  f  € 

j  \Ff(uri,(du)  <  ((/.  /»“/2  <  oo  (17) 

Proof.  Use  Lyapunov’s  inequality  (p.  191  of  Shiryayev,  1984)  and  (11)  to  get 
j^\Ff(uT^{du)sm\Ff(V)\‘  <  (E|F/(V)f)“^’  =  ((/,/))“/’“  <  «. 

The  proof  is  complete. 

Proof  of  Proposition  1.  The  functional  Ff{V)  =  Jq* e~*f{V{t))dt  =  F{f,V)  is. 
j<wtly  measureable  on  x  T.  Moreover,  for  ewJx  ^ed  /  €  Sj,  Ff(v)  €  L°(T)  (by 
Lemi^  5).  Thus,  (cp.  Rx>sinsld,  1986)  the  integral  in  (8)  is  well  defined,  ^eiefore 
$  has  the  integr^  representation  given  by  (8)  and  (cf.,  (3.2.2)  Samorodnitsky  and 
Tbqqu,  1990),  its  distribution  is  specified  by  (9).  It  remains  to  be  shown  that  $  is  a 
generalized  field.  Clearly,  $  is  linear.  By  Corollary  4.2  of  Walsh  (1986),  continuity  in 
probability  assures  that  the  field  has  values  in  the  dual  space  Sj. 

By  Theorem  2,  p.  254  of  Shiryayev  (1984),  and  the  linearity  of  it  is  enough  to 
show  that 

irj(t)  — »  0  in  as  j  -+  oo  ^  E\^(nj)\^  — » 0,  p  >  0.  (18) 

By  property  1.2.10,  p.  16  of  Samorodnitsky  and  Taqqu  (1990)  for  0  <  p  <  a 


=  (c(p.«))P  U 


where  c(p,  a)  is  a  constant  that  depends  on  a  and  p,  and  hence  does  not  affect  the 
convergence.  Thus  (18)  follows  immediately  firom  Lemma  5  and  the  proof  is  complete. 
Proof  of  Theorem  2.  By  Proposition  3.1  in  Szulga  (1992)  with  k  =  1,  it  is  clear 
that  as  n  — »  00,  $n(/)  converges  weakly  to  the  stable  integral  XFj  which  in  turn  is 
equal,  in  distribution,  to  f  €  S^}  (see  Proposition  1).  Use  the  Cramer- Wold 

device  (p.  49  of  Billingsley,  1^8)  and  the  linearity  of  $n  and  $  to  complete  the  proof. 

For  fh  €  Sj,  define  the  functional  •  •  • » I^)  as: 

^  ,  pj^)  _  I  F|^(Vi , . . . ,  Vi)  if  all  fc  arguments  are  distinct  ^ 

1 0  otherwise 


Since  the  probability  that  any  two  or  more  of  the  paths  of  the  Markov  processes  V  i , 
V^, . . . ,  Vn  are  identical  is  zero,  it  follows  that 


{4(Vi,...,Vi),  A  e sj}  2  A  e  sj}  (20) 


Consider  the  set  fig  =  {a;  :  Vi(w)  ^  ^  ^4(0;)}.  Clearly,  P(flo)  =  1-  Moreover, 

when  restricted  to  the  set  flo,  the  functionals  Ff^(Vi,. . . ,  I4)  and  Fy^(Vi, . . . ,  V4) 
are  identical.  Therefore,  for  the  next  two  results  in  this  section,  we  will  restrict  our 
atientimi  to  f2o>  end  denote  the  functionals  as  . . ,  V4). 

Proof  of  Proposition  3.  The  functional  Ff^iYi, . . . ,  V4)  is  jointly  measurable  on 

SjJ  X  T*.  Since  the  functions  /4(tti,...,U4)  =  /(ni)---/(nfc)  are  dense  in  5^,  it 
suffices  to  prove  the  Proporiticm  for  such  fig.  Note  that  in  this  case 


•  ••»«*)  =  Ff{ui) . . .  Ff{uig)  (21) 


Following  Kallenberg  and  Szulga  (1989),  we  set  C  to  be  a  Poisson  process  on  iZ\{0}  x  T, 
with  intensity  1/  x  /r,  so  that  (  is  constructed  the  jumps  of  the  process  X.  We 
denote  its  symmetrized  version  by  (.  Then  the  following  representation  of  the  integral 
holds 

■*■  27.  =  <  (^4)  (22) 


Here  L  is  the  operator  defined  on  the  space  of  functionals 


LFf^  =  LFf^(ui,. . . ,  tit;  zi, . . .  ,14)  =  zi . . .  X]gFf^{ui, . . .  ,U4), 


where  z,*  €  H  \  fO},  Uj  €  T,  i  =  1, . . . ,  fc  and  the  integral  on  the  right  of  (22)  exists  if 
Ff^  belongs  to  the  class  jC,  where 


£  =  {4  =  <*(^^4)^  <  “  »■'•}  ■ 

Cl) « *  • )  Ck  independent  copies  of  ^  and  note  that  by  Theorem  3.4  in  Kallenberg 
and  Szulga  (1989),  the  condition  ^^(LFf^  )^  <  00  a.8.  is  equivalent  to  Ci  •  •  •  Cki^^fk  ^ 
00  a.s.  By  the  independence  of  Cl)--*)Cifc)  (21)i  and  Lemma  2.2  in  Kallenberg  and 
Szulga  (1989),  this  is  equivaloit  to 


<  00  a.s. 


which  in  its  turn  follows  if  the  following  two  conditions  hold: 


0)  J  J  j  ^  v{dx)ii{dv)  <  00  a.s. 

(ii)  f  f  ({x^  A  1)(F?  V  1))  j/{dx)n{du)  <  00  a.s. 

JtJr\{0}  '  ^  ' 

Since  ^  l)v{dx)  <  00  for  Levy  measure  u  on  R\  {0}  and  (11)  holds,  (i)-(ii) 

are  satisfied.  Thxis,  the  multiple  integral  X^Fj^^  exists.  Linearity  of  in  /4  follows 
from. the  linearity  of  To  show  that  'i'  has  a  version  in  (<?|[)'  follow  the  proof  of 

Proposition  1,  keeping  in  mind  that  (cf.  (1.6)  of  Kallenberg  and  Szulga,  1989): 


=  x;  t!  (  *  )'  (v‘-“(xi...xt))'|  1^-”  (/-”■  (F;(ui)...F^(ui)))' 


where  /i(F/(u))2  =  Jj,  lF/(u)|2  fi(du)  =  (/,/)  <  c».  The  Proposition  is  proved. 
Proof  of  Theorem  4.  It  is  clear  from  Proposition  3.1  in  Szulga  (1992),  the  Crarner- 
Wold  device,  and  the  linearity  of  ^nifk)  that  the  finite-dimensional  distributions 

of  the  field  {^nifk)i  fk  ^  converge  weakly  to  those  of  ft  €  5^},  as 

n  — »  oo.  The  proof  now  follows  immediately  from  (20)  and  Proposition  3. 

Remark.  The  above  results  can  be  extended  in  two  directions.  The  first  one,  which  is 
qiiite  trivial,  is  to  introduce  the  parameter  ^  >  0,  to  substitute  the  Green’s  function  g 
iar  g  wherever  it  appeared,  and  to  change  the  exponent  c~*  to  in  the  defimtions 
(4)  and  (10).  We  will  denote  the  resulting  functionals  by  fj  and  the  spaces  of  the 

test  functions  by  5^’*,  s  5j.  It  is  clear  that  all  previous  results  hold  imder  these 

chanjm;  we  will  call  the  limiting  families  of  fields  {$^,0  >  0}  and  {^,0  >  0}. 

^e  above  results  also  hold  for  the  case  where  the  fields  are  indexed  by  mea¬ 
sures.  Let  us  briefly  explain  how  to  construct  Fjfj  =  Fyj(Jf,j, . .  .,Xi^)  when  7*  is  a 
symmetric  measure  on  ib  >  1,  such  that 


=  j  m(da)5^(a,x)/(x,y)7jfc(dx)7jt(dy) 


<  oo. 


(23) 


We  denote  the  class  of  such  measures  by  =  M.^ .  Theorem  2.1 

in  Adler  and  Epstein  (1987)  guarantees  existence  of  the  functional  F^^  for  each  7jt  G 
If  ^  is  absolutely  continuo\is  with  respect  to  Lebesgue  measure  with  density 
fk  then  .  Otherwise,  it  is  constructed  as  the  limit  of  path  integrals  of 

the  fonn  g  :=  . Vt(ti))*i  Here,  ifidx)  = 

=  Thus.^f 

is  a  smoothed  version  of 


4.  POINT-INDEXED  STABLE  RANDOM  FIELDS 


In  the  previous  section  we  have  shown  the  construction  of  the  stable  random 
fields  indexed  by  functions  /  €  Sjj.  We  will  now  show  how  the  same  construction 
allows  us  to  obtain  stable  fields  on  To  do  so,  we  would  like  to  apply  Theorem  2  to 

the  case  of  measure  7i(dz)  =  Sx{z)dz,  where  x  €  R^.  Of  course,  this  is  possible  only 
when  7x  €  A<^,  i.e.,  when  (23)  holds.  This  is  true,  for  example,  for  Brownian  motion 
on  R^.  The  corresponding  functional  Fy^{V)  is  then  the  exponentially  weighted  local 
time  Lx(V)  of  the  process  V  at  point  x. 

Denote  by  x  €  R^,  the  stable  random  field  obtained  by  applying  the  measure- 
variant  of  Theorem  2  to  the  sums  of  the  functional  LxiY).  Then  x  €  R^,  h*is  the 
integral  representation 

$(x)  =  ix  =  J  Lx{u)X{du). 

So,  whenever  the  local  time  of  the  Markov  process  exists,  we  obtain  the  point-indexed 
stable  random  field  {^(x),  x  G  R!^},  x  G  with  finite-dimensional  distributions 
given  by 

Eexp{»  ti^xi }  =  exp  I -E  U^XiiV) 

t=l  I  li=il 

Clearly,  when  a  point-indexed  field  {^(*),  x  €  exists,  we  can  always  create 
a  indexed  version  of  it  by  setting  $(/)  =  $(x)/(x)dx. 


In  a  similar  manner,  for  those  processes  for  which  measure  = 

~  X2)dx\dx2  is  of  finite  norm,  one  can  use  exponentially  weighted  intersec¬ 
tion  local  times  I'xCVi,  Vjj)  to  define  the  point-indexed  random  fields  x  €  /?^} 

as  limits  of  sums  V-,).  By  Theorem  4 

®(®)  =  J^^-^'x(«i,U2)Jf(dui)X(du2). 


S.  RENORMALIZABLE  STABLE  FIELDS 

The  fxmoept  of  zenormalizability  of  families  of  generalized  random  fields  was  in¬ 
troduced  in  Adler  and  Epston  (1987^  as  an  extennon  of  the  property  of  sdf-similarity. 
There  the  authors  considered  conditions  for  existence  of  this  property  for  some  Gaus¬ 
sian  families  and  their  multiple  integrals.  Here  we  modify  th^  definition  of  renor- 
malizability  to  include  rescaling  of  sn  additi<mal  parameter  and  study  this  property 
for  the  famihes  of  stable  fields  that  are  obtained  as  limits  of  functionals  of  the  paths 
of  self-similar  Markov  processes. 

In  the  following,  we  will  assvune  that  the  imtial  measure  of  the  Markov  process 
is 

(24) 

where  B  C  is  a  fixed  set  of  finite  Lebesgue  measure  \B\. 

Definition.  A  process  V  is  said  to  be  self-similar  with  index  /3,  if  for  any  ij  >  0,  V 
and  the  process 

(25) 

are  identical  in  distribution,  where  V  is  a  process  (on  another  probability  space)  which 
has  the  same  transition  probabilities  as  V  but  initial  measure  m(dx)  =  j^l^ff(x)dx. 

Note  that  the  initial  measure  has  to  be  changecT  and  the  set  B  has  to  be  scaled  in 
order  to  preserve  the  distribution  of  the  starting  points  after  scaling  V  by  the  factor 
t}.  For  any  function  /(x)  on  and  t  >  0,  define 

^/(x)  :=  i7~V(i?-^x).  (26) 

Lemma  6.  Let  V  be  a  Markov  process  which  is  self-similar  with  index  P,  and  let 
q  €  Sj’*,  jb  >  1.  Then 


^(Vi . . %),r=k?  (2?) 

Proof.  Using  the  self-similarity  property  of  V,  and  then  malcing  the  change  of  vari¬ 
ables  tj^ti  5=  Si,  we  get 

=  /j, . ...it* 


O 


with  T  =  k$.  This  proves  the  lemma. 

Let  us  denote  by  the  SaS  Levy  random  measure  X(.),  on  (T,  T ,  /i*”),  with 

control  measure  /i”*.  The  index  m  indicates  that  the  process  V  starts  according  to 
the  initial  measure  m.  FVirther,  let  denote  the  multiple  integral  of  order  k,  of 

the  ftmction  /  with  respect  to  When  fc  =  1,  we  set  =  X”*/.  We  define 

the  following  two  families  of  stable  fields: 


Definition.  The  fainily  of  random  £elda  d  >  0}  is  renormalizable  with  renor¬ 

malizing  parameters  (r,  p)  if,  for  every  17  >  0, 

V 

where 

»7  \V  /  V 

If  the  are  measure  indexed,  set 

>7  \»7  /  »7 

Theorem  7.  Let  (V(f),m)  be  a  self-similar  Markov  process  with  index  and  let 
g^,  ff  >  0,  be  the  corresponding  Green  function.  The  family  of  stable  random  Gelds 
>  0}  is  renormalizable  with  parameters 
Proof.  Using  (9)  and  Lemma  6  with  r  —  fi,  one  obtains; 

Eexp  lit**’''’*  (^,)  }  =  exp  |-rE  |i^''(V)|°} 

=  exp{-|«rE|j^(V)|“} 

=  Eexp {*<$*’’"•(<7)}  ,  teR 

and  the  proof  is  complete. 

Our  next  result  establishes  the  renonnalizability  of  the  family  ^  >  0} 

of  mtiltiple  integrals  of  the  fields  0  >  0}. 

Theorem  8.  Let  fc  >  1.  The  family  0  >  0}  is  renormalizable  with  parameters 

Proof.  Let  “=j^”  denote  weak  convergence.  By  Theorem  4, 

J2  -  -  »sn-.oo, 

1<*1<— <**<n 

X)  y-i,  ■  ■  ■  J'm.  .  ■  •  • .  Vi, )  =*•  ‘fgi)  as  n  -  00. 

By  Proposition  3,  ’^^•"*(9*)  =  qt)  =  and  by 

1)  9* 

Lemme  6,  Jij( Vi, . . . ,  Vj)  £  fff ' (Vj, . . . , Vi)  Thus.  we  have. 

«?  ® 


and  the  theorem  is  proved. 


6.  SIMULATION  OF  STABLE  RANDOM  FIELDS 

Here  we  provide  graphs  of  some  J^^-indexed  stable  fields  for  d  =  1  Jind  2.  The 
aunulations  were  based  on  results  of  Theorem  2  and  discussion  in  Section  4,  which 
allow  one  to  btdld  the  fields  firom  the  stuns  of  functionals  of  the  Markov  processes 
V.  We  would  like  to  explain  briefly  how  these  simulations  were  done. 

We  simulate  symmetric  random  walks  St  on  the  integer  lattice  Z**,  d  =  1, 2,  i.e., 

5o  =  0,  +  ...  +  (*,  k  =  1,2,-.., 

where,  the  (i  are  i.i.d.  Rademacher  random  variables.  Their  initial  measure  is  m(dx)  = 

l£(x)dz,  where  S  C  is  a  fixed  set  of  finite  Lebesgue  measure  \B\  to  be  specified 
later.  Following  a  discussion  in  Oynldn  (1988),  use  the  functional  of  the  random  walk 

tnT 

i=l 

to  approximate  the  value  of  the  exponentially  weighted  local  time  £/  of  the  Brownian 
motion  (which  exists  for  d  =  1).  This  is  a  go^  choice  for  several  reasons:  As 
m  -*  oo,  the  finite-dimensional  distributions  of  the  scaled  random  walk  ^  ^ 

{0, 1/m,  2/m, . . .}}  on  the  lattice  Z^|^/m  converge  to  those  of  the  Brownian  motion 
on  We  would  lilcw  to  have  a  functional  of  S/ y/rri  which  approximates  the  functional 

of  Brownian  motion 

FfiV)  =  jf  °°  (35) 

Since  the  value  e"^  is  less  than  10“®,  for  computational  purposes  we  restrict  the 
interval  of  integration  in  (35)  to  (0,  T),  with  T  =  20.’Dynkin  (1988)  gives  the  approx¬ 
imation  as 

Tm 

i=l 

In  order  to  obtain  point-indexed  fields,  we  replace  the  function  /  by  the  Dirac  delta- 
function  at  1.  Finally,  we  use  smns 


<=1 

to  approximate  values  of  the  stable  field  $(/). 

For  the  case  d  =  1,  we  have  chosen  m  =  1000,  n  =  2000,  and  B  =  (—1000, 1000]. 
We  circularized  the  interval,  i.e.,  when  the  random  walk  wanders  off  the  left  hand- 
side  of  [—1000, 1000],  it  immediately  reappears  at  the  right  hand-side  and  vice-versa. 
To  fix  the  value  of  m  we  first  obtained  graphs  of  (34)  for  one  fixed  realization  of 
the  random  vralk  (but  different  values  of  m)  and  concluded  that  the  graph  does  not 
change  much  for  m  >  1000.  To  choose  the  v^ue  of  n  we  fixed  a  =  1.95  and  m  =  1000 
and  generated  independently  n  =  1500  values  of  the  stable  process  at  a  particular 
pcnnt.  This  histogram  was  very  similar  to  the  histograms  of  1500  values  from  a  SaS 
distribution  with  a  =  1.95  and  a  —  0.014678  generated  by  the  software  package  Splus. 
We  also  estimated  the  values  of  the  parameters  a  and  a  using  a  method  based  on 
McCullough  (1986).  The  value  of  a  was  close  to  1.95  and  a  =  0.014678. 

For  d  =  2  we  took  m  =  1500,  n  =  4000  and  B  =  [0, 100)  x  [0, 100).  As  for  d  =  1, 
we  drcularized  the  region. 


in 


To  generate  the  stable  random  variates  that  satisfy  (2)  we  followed  the 
algorithm  in  Bratley,  Fox  and  Schrage  (1987),  which  generates  SaS  random  variates 

with  scaling  factor  a  =  1,  and  then  multiplied  the  values  by  The  value  of  a  is 

chosen  before  simulating  the  stable  and  Rademacher  random  variates. 

The  results  of  our  simtilations  are  presented  in  Figures  6.1-6.6.  Figures  6. 1-6.3 
present  SaS  processes  ♦(x),  x  =  //m,  I  =  -1000,  -999, . . . ,  1000  for  three  values  of 
the  parameter  a  =  1.1, 1.5  and  1.95.  The  points  1  =  0  and  2001  on  the  graphs  in  the 
Figures  correspond  to  the  points  1  =  -1000  and  x  =  1000,  respectively.  The  graphs 
of  the  stable  fields  ^(x),  x  €  [0, 100)  x  [0,100)  for  a  =  1.1, 1.5  and  1.95  are  given  in 
Figures  6.4,  6.5,  and  6.6.  Elam  of  tnese  figures  contains  three  graphs.  The  graph  at 
the  top  is  the  graph  of  the  stable  random  fidid.  The  graph  at  the  bottom  left  is  of 
the  positive  values  of  tiie  fidd  plotted  separatdy;  and  the  gri^h  at  the  bottom  right 
is^ of  the  ne^tive  values  of  the  field.  It  is  dear  fiom  the  grap^  that  when  a  is  close 
to  1,  the  fidd  is  &irly  flat;  as  a  increases  the  fidds  become  more  peaked. 
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Fig.  6.4.  Stabk  Rudom  Field,  a  ss  1.1  :  (»)  field,  (b)  podtive  veluee  of  the  field,  (c)  negative  value*  of  the  field 


Stable  Random  Field,  a  as  1.S  :  (a)  Held,  (b)  pocitive  values  of  the  field,  (c)  negative  values  of  the  field 


SUbk  Random  Field,  <k  s  1.05  :  (a)  Ikld,  (b)  podtive  valuce  of  the  field,  (c)  negative  values  of  the  field 
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